Poincare gauge theory of gravity: Friedman cosmology 
with even and odd parity modes. Analytic part 



Peter Baekler* 
Department of Media, Fachhochschule Diisseldorf, 
University of Applied Sciences, 404^4 Diisseldorf, Germany 

Friedrich W. Hehl^ 

Institute for Theoretical Physics, University of Cologne, 50923 Koln, Germany, 

and 

Department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA 

James M. Nested 

Department of Physics & Institute of Astronomy & Center of Mathematics and Theoretical Physics 
National Central University, Chungli 320, Taiwan, ROC 

(Dated: 06 November 2010, file PGCosmology/PGcosmologyOne40X.tex) 

We propose a cosmological model in the Iramework of the Poincare gauge theory of gravity (PG). 
The gravitational Lagrangian is quadratic in curvature and torsion. In our specific model, the 
Lagrangian contains (i) the curvature scalar R and the curvature pseudo-scalar X linearly and 
quadratically (including an RX term) and (ii) pieces quadratic in the torsion vector V and the 
torsion axial vector A (including a VA term). We show generally that in quadratic PG models 
we have nearly the same number of parity conserving terms ('world') and of parity violating terms 
('shadow world'). This offers new perspectives in cosmology for the coupling of gravity to matter 
and antimatter. Our specific model generalizes the fairly realistic 'torsion cosmologies' of Shie- 
Nester-Yo (2008) and Chen et al. (2009). With a Friedman type ansatz for an orthonormal coframe 
and a Lorentz connection, we derive the two field equations of PG in an explicit form and discuss 
their general structure in detail. In particular, the second field equation can be reduced to first 
order ordinary differential equations for the curvature pieces R(t) and X(t). Including these along 
with certain relations obtained from the first field equation and curvature definitions, we present 
a first order system of equations suitable for numerical evaluation. This is deferred to the second, 
numerical part of this paper. 

PACS numbers: 04.50.Kd, 11.15.-q, 11.30.Er, 98.80.Jk 



I. INTRODUCTION 

In order to accommodate the local Poincare group in 
spacetime, Sciama [1] and Kibble [2] had to extend the 
Riemannian spacetime of general relativity (GR) to a 
Riemann-Cartan spacetime with non-vanishing torsion 
T a (for the notation see the end of the Introduction). 
Thereby the orthonormal coframe i9 a and the Lorentz 
connection F Q ^ = — T^ a became independent gauge po- 
tentials of weak and strong gravity, respectively. The cor- 
responding gauge field strengths are torsion T a = D§ a 
and curvature R a/3 ~ DT a P , as spelled out in Sec. II. 
There we also display the irreducible decompositions of 
T a and R af3 . 

If one allows in a Yang-Mills manner for a gravitational 
Lagrangian V that is quadratic in torsion and curvature, 
we speak of a Poincare gauge theory of gravity (PG) [3- 
10]. In Sec. Ill A we introduce the gravitational excita- 
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tions H a = -dV/dT a and H af} = -dV/dR a ^ and reca- 
pitulate the general form (21) and (22) of the two field 
equations of gravity. 

Then, in Sec. IIIB, we turn to the conventional parity 
conserving quadratic Lagrangian V+, which includes the 
somewhat degenerate Einstein-Cartan Lagrangian Vec- 
Because of the existence of the Euler 4-form of the cur- 
vature, we can show that one curvature square piece is 
trivial. In Sees. Ill C and HID, we review parity vio- 
lating admixtures to the EC-Lagrangian that have been 
formulated in the past by different groups, stressing in 
Sec. Ill E the importance of the corresponding cosmolog- 
ical models of Shie-Nester-Yo [ ] and Chen et al. [ ]. 

Having in this way the PG at our disposal, we open 
for it in Sec. IV a new 'window' to a 'shadow world': In 
Sec. IV A, we show in a systematic way that, besides the 
parity conserving Lagrangian V+, there exists an equally 
important Lagrangian V- the pieces of which are all par- 
ity violating. Accordingly, for PG we propose the grav- 
itational Lagrangian V± = V+ + V— . An equivalent La- 
grangian has already been discussed earlier by Obukhov 
et al. [14]. 

Because of the complexity of this general Lagrangian, 
we select for further study in Sec. IV B in Eq.(64) the sim- 
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pier 9-parameter Lagrangian Vbhn, which should carry 
the characteristic features of parity conserving and par- 
ity violating effects. In Sec. IV C a novel method is pro- 
posed for diagonalizing the quadratic pieces in Vbhn ■ No 
linearization is involved and the output consists of exact 
analytic results. Besides the Einstein mode 2 + , we find 
for the torsion modes spin and parity , and for the 
curvature modes Q ± . 

We calculate the gravitational excitations of Vbhn 
(Sec. IV D) and display in Sec. IV E the corresponding 
field equations explicitly. They turn out to be first order 
partial differential equations in torsion and curvature, re- 
spectively. 

We continue by looking closer into the structure of 
Vbhn and its field equations. In Sec. IV F the Nieh-Yan 
identity is used to show that the coupling constants of 
Vbhn occur only in certain linear combinations in the 
field equations. 

Eventually, in Sec. V, we turn to a cosmological model. 
In Sec. VA the coframe and the torsion are assumed 
to be homogeneous and isotropic in accordance with 
a Friedman-Lcmaitre-Robertson- Walker (FLRW) type 
model and in Sec. VB the corresponding irreducible 
pieces of the curvature are calculated. We define a spin- 
less perfect fluid in Sec. VC and find then, in Sec. VD, 
the field equations of gravity for this cosmological model. 
The first field equation yields equations for the density 
p(t) and for the pressure p(t) of the perfect fluid. These 
equations are subsequently manipulated in order to bring 
them into a more transparent form. The second field 
equation has also two independent components, namely 
first order ordinary differential equations for R(t) and 
X(t). We uncouple them and bring them in the very 
compact form (176) and (177) by introducing certain 'fre- 
quencies' Wo, CJl, UJ2, W3. 

Now we are able to evaluate our exact results by nu- 
merical methods. This will be done in follow up work. 



Notation 

Our notation is as follows (see [15, 16]): We use the 
formalism of exterior differential forms. We denote the 
frame by e a , with the anholonomic or frame indices 
a, /?, • • • = 0, 1, 2, 3. Decomposed with respect to a natu- 
ral frame di, we have e a = e l a di, where ■ ■ ■ — 0, 1, 2, 3 
are holonomic or coordinate indices. The frame e a is the 
vector basis of the tangent space at each point of the 4D 
spacetime manifold. The symbol j denotes the interior 
and A the exterior product. The coframe d@ = e^dx^ is 
dual to the frame, i.e., e Q ji9^ = 5^, The * denotes the 
Hodge star operator that acts on the quantities on its 
right, as, for instance, in *(E a A ¥). If d aP := iT A ¥ , 
etc., then we can introduce the eta-basis by 77 := *1, 
i] a := e a jr) = *iT, rf^ := = *-d af \ etc.. Paren- 

theses surrounding indices (a/3) := (a/3 + /3a) /2 denote 
symmetrization and brackets [a/3] := (a/3 — /3a) /2 anti- 
symmetrization. 



The coframe ~d a and the 77-system are related by 

tT a mi = 5° v fi - s<$ 777 , 

-d a A n^s = £5 ?7 ( 3 7 + 5% risp + 8% r) yS , (1) 

Differentiating the 77's, we find in a metric-affine space 
(for a definition see Sec. II B) the relations 

Drj a = 2Q A rj a + A rj afJi , 

Driufj = —2Q A rj a p + A rj af3fl , 

D-q a01 = -2Q A r; Q( 3 7 + T M A 77^-^ , 

Dria^s = ~2Q A T7 Q/37<5 , (2) 

where Q is the Weyl covector and the torsion. 

We use the abbreviations: GR = general relativity the- 
ory (Riemann spacetime), EC = Einstein- Cartan theory 
of gravity (Riemann-Cartan spacetime with torsion and 
curvature, gauge Lagrangian linear in curvature), PG 
= Poincare gauge theory of gravity (Riemann-Cartan 
spacetime, gauge Lagrangian arbitrary function of tor- 
sion and curvature; often a quadratic function) , MAG = 
metric-affine theory of gravity (metric-affine spacetime, 
gauge Lagrangian arbitrary function of in torsion, non- 
metricity, and curvature; often a quadratic function). 

II. SPACETIME OF GAUGE THEORY OF 
GRAVITY 

A. The coframe and weak gravity 

One of the fundamental structures in a gauge theory 
of gravity is the coframe field t9 q — and all quantities are 
referred to it. It is represented by four linearly inde- 
pendent 1-forms 7? a , with a = 0,1,2,3. They can be 
decomposed with respect to a natural coframe dx 1 ac- 
cording to t9 q = ei a dx % . Here the Ci a are the coordinate 
components of also called tetrad components. 

Since a Riemannian metric g with Lorentz signature is 
assumed to exist, the coframe can always chosen to be 
orthonormal according to 

g = gij dx i <g) dx j = g ajj tT ® -d 13 , (3) 

with g a p — diag(— 1, 1, 1, 1). In this way the metric is 
absorbed by the coframe and has no longer independent 
physical degrees of freedom. 

This choice is convenient and will be kept through- 
out this article. However, sometimes one may want to 
choose arbitrary coframes: then the metric emerges ex- 
plicitly again. Also for that reason, the metric and the 
coframe, besides the linear connection (see Sec. II B), are 
treated in the variational principle as independent gauge 
field variables. However, it eventually turns out that the 
field equations resulting from the variation of the met- 
ric and of the coframe are equivalent. In other words, 
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the ortho normal "gauge" of the coframe, which we use in 
this article, doesn't restrict the generality of our consid- 
erations. 

We call d a the potential of weak gravity of the Newton- 
Einstein type. It couples to matter via the Einstein grav- 
itational constant k, which has the dimension of a re- 
ciprocal force. Basically, d a represents four gauge boson 
fields — each of helicity 1 — that conspire, at least in linear 
approximation, to build up the massless spin 2 graviton 
modes of general relativity (GR), provided the appropri- 
ate Hilbert Lagrangian is chosen. 



B. The linear connection Y a and the hypothesis of 
strong gravity 

In classical gauge theories of gravity the linear connec- 
tion of spacetime T Q ^, besides the coframe $ a , is assumed 
to exist as an independent field variable. This is a phys- 
ical hypothesis that has eventually to be checked by ex- 
periment. We call T a " the potential of strong gravity of 
Yang-Mills type. It is represented by 4 x 4 = 16 bosonic 
one-forms fields that can be decomposed according to 
T a P — Tia^dx 1 . They couple to matter in a Yang-Mills 
like fashion via a hypothetical coupling constant g of the 
dimension (action) . 

A differential manifold equipped with a metric g and 
a linear connection T a ^ is called a metric-afnne space. 
If the linear connection is unconstrained, then r ct ' 9 has 
values in the Lie algebra of the general linear group 
GL(4,R). The gauge theory with independent metric 
and independent connection is called metric- affine (gauge 
theory of) gravity (MAG). Since in MAG the connection 
components Ti a ^ carry three indices, the strong gravity 
potential can provide additional strong gravity modes of 
up to spin 3 (see [17]). 



C. Beyond general relativity: relaxing the torsion 

Let us start from GR and assume the connection to 
be Riemannian, namely T a ^; we will always indicate 
the Riemannian nature of a quantity by a tilde. Then 
T a @ doesn't provide a mode that is independent of the 
coframe. With a suitable Lagrangian, this corresponds 
to GR. 

If the connection is metric-compatible but non- 
Riemannian, it carries an independent piece that has val- 
ues in the Lie- algebra of the Lorentz group SO (1,3). In 
orthonormal frames, we have for the Lorentz- (or spin-) 
connection the relation L Q/3 = — T@ a . It represents 6 
bosonic one-form fields for strong gravity. Its maxi- 
mum spin is 2. The independent fields, coframe "d a and 
Lorentz-connection T Q ^, represent the gauge potentials 
of the Poincare group. The corresponding gauge field 
theory is called the Poincare gauge theory of gravity 
(Poincare gravity or PG). 



In PG, the contortion 

K a/3 •= f Q ^ — r a/3 = —K^ a (4) 

measures the difference between the Riemann and the 
Riemann-Cartan (RC) geometry and, as a difference be- 
tween two connections, it constitutes a tensor. Alterna- 
tively, the deviation from the Riemannian geometry of 
GR can be described by the torsion 

T a :=D 0° = M a + T a A & = ^T ij a dx i A dx j . (5) 
r 

Here D is the exterior covariant derivative with respect 
to the connection T a @. The newly emerging Lorentz- 
connection modes reflect themselves also in those of the 
torsion, since the second term in the torsion, Tp a A d^ 3 , 
depends on T a ° . It can be shown (see [15]) that torsion 
and contortion are related by 

T a = A Kp a , K a0 = e [a jT p] - X - (e QJ e^T 7 )iF. (6) 

On the level of the gauge potentials, we have then in 
PG the frame i? Q and the Lorentz-connection r a/3 = 
— r' 3 ". On the level of the newly introduced torsion, 
the translation gauge field strength, we can execute an 
irreducible decomposition in order to learn more about 
its structure. It can be decomposed according to 24 = 
16 © 4 © 4 into three pieces: into a second rank tensor 
piece (tentor), Wt q , into a vector piece, namely, in com- 
ponents, the trace of the torsion (trator), 

(2) T a . = VAir with V:=e (3 jT' 3 , (7) 
o 

and into an axial vector piece (axitor) , which corresponds 
in components to the totally antisymmetric piece of the 
torsion (the star denotes the Hodge operator): 

(3) T Q = \*{AA $ a ) with A := * (?9 Q A T a ) . (8) 
We have then the irreducible decomposition 

tentor trator axitor 

The tensor piece can carry at most spin 2 modes, whereas 
the vector and the axial vector pieces are good for at most 
spin 1 modes. 



D. Keeping the nonmetricity to zero 

In MAG, besides the torsion T a , we have a nonvanish- 
ing nonmetricity 

r 

Q a p '•= — D g a /3 = Qp a = -dg a p + 2r (ct(3 ) . (10) 
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Therefore, the nonmetricity Q a p of the metric-affine ge- 
ometry is a measure for the difference between the linear 
connection Y a @ and the Lorentz-connection. In an or- 
thonormal coframe i9 Q the metric referred to the coframe 
g a /3 is a constant and the modes of the symmetric 1 
rVa/3) — ^i(aj3)dx l are passed through, according to (10), 
to the nonmetricity Q a p = Qiapdx 1 , with the compo- 
nents Qi a fi = Qifia- There emerge, besides the 6 of the 
Lorentz-connection, 10 more bosonic one-form fields Q a p 
clearly encompassing strong gravity contributions of spin 
0,1,2, and 3. 

In MAG, we have the gravitational potentials $ Q and 
T a P . The gauge field strength attached to the coframe i} a 
is the torsion (5), that attached to the linear connection 
the curvature 2- form 



R, 



= dT p — r 



, 7 a r/ . 



(ii) 



We can raise the index a and can decompose the cur- 
vature into the antisymmetric 'rotational' piece and the 
symmetric 'strain' piece according to 

R a P = W al) + Z a/3 (12) 
with W afi := R [af3] , Z a(i := . 

Even though we succeeded to relate a quadratic MAG- 
Lagrangian to a consistent classical field theory of mass- 
less spin 3 fields via the tracefree part of the nonmetricity 
[17], we will restrict ourselves in this article to vanishing 
nonmetricity, Q a p = 0; consequently the strain curva- 
ture vanishes, too: Z a/3 — 0. From a phenomenolog- 
ical point of view the overwhelming importance of the 
(rigid) Poincare group in special relativity directs our at- 
tention primarily to the gauge theory of the local Poincare 
group, namely PG. Accordingly, from now on Z a @ = 
and R afs = W afs . 



E. Irreducible decomposition of the rotational 
curvature 

For the physical interpretation it is significant to un- 
derstand the different pieces of the curvature. The rota- 
tional curvature decomposes irreducibly into six pieces, 
see [15, 18-21], according to 



R, 



a/3 



(1) R a ,3 + ^Ra/3 + i3) R a p 
paircom 9 

ricanti 6 



weyl 10 

{4) R a p 



pscalar 1 

{6) R a p . 



(13) 



ricsymf 9 



scalar 1 



The number of independent components is specified sub- 
sequent to the (computer) name of the corresponding ir- 
reducible piece. Pseudoscalar and scalar qualify as linear 



1 Not to be confused with the symmetric part of the connection 

r( a ^) 7 = e( ct jr /3 ) 7 . 



Lagrangians. We take from the literature (remember that 
$ a p =iA^ and r] a p = *$ a p) 



&R, 



a0 



1 

12 



(14) 



X:=e a jX a , X a :— * A dp) ; 

(6) R*p = ~R$ a p, (15) 

R := e a jR a , R°> := epjR a0 . 

We recognize that the scalar R and the pseudoscalar 
X play a preferred role. Note that X is purely post- 
Riemannian, that is, X = 0. In components we have 

X = VaPl5 R^ l3 ^/i\ and R = R f3a al3 , (16) 

with the decomposition R a/3 = R lS al3 -d~ <5 /2. 



III. POINCARE GAUGE THEORY OF 
GRAVITY (PG) 

A. Lagrangian and field equations 

In PG, we have the gravitational potentials -d a and 
r Q/3 = — r^ Q . The corresponding gauge field strengths 
are the torsion T a and the rotational ('Lorentz') curva- 
ture 

Ra p ._ dT ap _ A p 7 /3 = _ R p a ( 17 ) 

We assume a first order Lagrangian consisting of a gauge 
and a minimally coupled matter part, 



L = V (g a0 , d a , T a , R al3 ) + L mat (g a p , iT , * , D * 



(18) 

with the matter field(s) '3/ . Then we can define the trans- 
lation and the Lorentz excitations, respectively, 

dV dV 
~ ~ qj^ ) H a p;=- 7r ^ = -Hp a , (19) 



dR a P 



and the canonical matter currents of energy-momentum 
and spin (angular momentum) according to 



-rp a ; (20) 



in the case of a minimally coupled matter Lagrangian, as 
in (18), the variational derivatives degenerate to partial 
derivatives. 

The action principle yields the field equations [3] 

DH a — E a = E Q (first), (21) 
DH aj3 - E aP = T af3 (second) , (22) 

with the gauge currents of energy-momentum and spin 

E a := e a jV + (e^T* 3 ) A H + (e a jR^) A H M , (23) 



E, 



a/3 



-i?r„ A H, 



P] 



(24) 
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If the gauge Lagrangian V is prescribed explicitly, we 
can compute first the excitations H a , H a p by partial dif- 
ferentiation of V and subsequently the gauge currents 
E a ,E a p by substitution; these quantities are then in- 
serted into the two field equations (21), (22). As noted 
already above, the field equation resulting from a varia- 
tion of the metric g a p is equivalent to (21), provided (22) 
is fulfilled. 

The matter currents on the right-hand-side of the field 
equations (21) and (22) can be understood as those of a 
spin fluid (see [22-24]). An approximate representation 
of such a spin fluid can be specified as follows: If the 
fluid moves with the velocity u — u a e a , that is, with the 
flow 3-form U := ujv = u a rj a , and transports an energy- 
momentum density p a and a spin density s Q/ g = — s^ Q , 
then a convective Weyssenhoff ansatz for the matter cur- 
rents reads 



p a U and r, 



a/3 



S a/ 3 



u. 



(25) 



B. 



Quadratic Yang-Mills type Lagrangian with 
even parity terms 



A quadratic Lagrangian of PG of the Yang-Mills type 
has the general structure (Ao is the cosmological con- 
stant) 

V ~ — (curv + torsion 2 + An) + - curv 2 . (26) 

K Q 

Since the coframe r d a and the Lorentz connection T a/3 
are independent variables, such a first order Lagrangian 
yields second order field equations; higher derivatives do 
not emerge. 

The simplest non-trivial Lagrangian corresponds to the 
first term on the right-hand-side of (26). It is of the 
Hilbert type, i.e., linear in the curvature, namely the so- 
called Einstein-Cartan Lagrangian, see (15), 

Vec = ^V a p A R a ' 3 = A ^R a ? = ±-*R . (27) 

Ait\, 

The corresponding two field excitations, H a = and 
H a p = — r) a p/(2K), if substituted into (21) and (22), yield 
the field equations [1, 2, 25]: 



-TfaP-y AT 7 = KT a p . 



(28) 
(29) 



The viable Einstein-Cartan(-Sciama-Kibble) theory 
(EC), as compared to GR, supplies an additional spin- 
contact interaction of weak gravitational origin, since 
only Einstein's gravitational constant enters (28), (29). 
The Lorentz connection T"' 9 cannot propagate and thus 
EC represents 2 a degenerate PG. In order to enable T a " 



to propagate, we have to use additionally at least the 
quadratic curvature piece in (26); for discussions on the 
physical relevance of torsion one should also compare 
Shapiro [ ] and Ni [29]. 

Esser [30], see also [31], constructed the most general 
quadratic Lagrangian with even (+) parity pieces (for 
this notion see Sec. IIIC). For a RC-spacetime it reads: 



2 k 



{-a R a P A r, aP - 2X r, 

+T a A £ aj^T a ) 
1=1 



Z Q i=i 



(30) 



without restricting the generality of our considerations, 
we can choose g > 0. In a RC-space, Esser found the 
additional term 



— R a/3 A * 
2g 



w 7 <& a A (e 7 J (5) i? 7 /3 ) 



(31) 



However, this term can be transformed successively into 
a pure 105 term: Let us first consider invariants of the 
form R afi A * [i? Q A (e 7 _i ^JF^)] with A e {1---6}. 
For A = 5 we find 



R a P A * 



.< B >JF, 



d a A [e 1 A 
= &R a0 A * [tf a A (e 7 J (5) i?T^) 
= (5) i? Q/3 A * (5) R af3 . 



(32) 



Thus we can absorb the wj term into the w$ term. Con- 
sequently, without restricting the generality of our con- 
siderations, we can put wj = 0. 

In our Lagrangian V+ in (30) not all the constants 
are independent. In a Riemannian as well as in a RC- 
spacetime, the integrand of the topological Euler 4-form 



B 



Hill') - -^Ra 13 A R^ a - -r/fSas^ 



R aP A R' (S 



(33) 



= dC RR <,) 
is exact, with 

Crrm = \r\ a & \ (iz^Ar^+ir^AiyAr/) . (34) 

The dual is here taken with respect to the frame indices 
a,/3 of the curvature 2-form R a ^, it has to be carefully 
distinguished from the Hodge dual. Then, together with 
the Bach-Lanczos identity [15], Eq.(A.3.7), 



flH7 Aj R(*)l« 7 _ lg°0 Rl»> A r{*) =0 , 



(35) 



2 Recently, an EC-model with fermionic matter and its application 



to the early universe has been discussed by Ribas & Kremer [2 
and Dolan [27]. 
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one can show that only five of the six ui/'s are linearly 
independent. 

The excitations can now be calculated by differentia- 
tion: 

fla = -- E aj*«T a , (36) 

i?a/3 = %a/3 + - E «)/ * (/) ^ . (37) 

2k g /=1 

Because of (15), the last equation can be slightly rewrit- 
ten as 

H aP =(^-^R)va0+ 1 Em^Ra0. (38) 

\2k 12 g J H g I=1 

There have been numerous investigations into the proper- 
ties of the Lagrangian (30). In linear approximation, on 
a flat Minkowskian background, Eq. (30) encompasses, 
besides the weak gravity modes of the coframe, propa- 
gating strong gravity modes of the Lorentz-connection 
with spin 2 ± , l*, and O*, as shown by Hayashi & Shi- 
rafuji [ ], by Sezgin & van Nieuwenhuizen [32], and by 
Kuhfuss & Nitsch [33]. For a model with quadratic cur- 
vature Lagrangian in which only the Lorentz-connection 
is dynamic, compare Cho et al. [34]. 

A good dynamic mode transports positive energy at 
speed < c. At most three modes can be simultaneously 
dynamic; all the cases were tabulated; many combina- 
tions are satisfactory to linear order. The Hamiltonian 
analysis, as shown by Blagojevic and Nikolic [35, 36], 
revealed the related constraints. In more detailed in- 
vestigations [37-40] it was concluded that effects due to 
nonlinearities could be expected to render all of these 
cases physically unacceptable, with the exception of two 
"scalar" connection modes with spin + and spin 0~. 

Before we come back to the mode analysis in Sec. IV C, 
we want to extend the gravitational Lagrangian such that 
also odd parity pieces are included. 

C. Even parity and odd parity Lagrangians, 
twisted and untwisted forms 

Let us study the spatial reflection or parity transfor- 
mation; the sign of the time axis will be kept fixed. A 
(pure) scalar field remains invariant under the par- 
ity transformation or if we transform a right-handed co- 
ordinate system x % into a left-handed one x % : <&(x') = 
<j>(x). In contrast, a twisted scalar field $(z) (also called 
a pseudoscalar field) changes its sign under those cir- 
cumstances, that is, the sign of the determinant J := 
det 1 1 dx 1 /dx 1 \ | of the Jacobian transformation matrix en- 
ters its transformation law: <t(x l ) — (sign J)&(x). The 
analogous behavior characterizes the relation between 
twisted and untwisted forms; for a mathematical discus- 
sion compare Frankel [41]. 

A Lagrangian 4-form L has to be a twisted 4-form in 
order to make its action W := J n La pure scalar. These 



Lagrangians are also called even parity Lagrangians. 
However, in physics we know since the discovery of parity 
violation in the weak interaction in 1956, see Sozzi [42] 
for a review, that also odd parity would-be Lagrangians 
can occur; they have to be multiplied by pseudoscalar 
coupling constants in order to transform them to decent 
(twisted) Lagrangians, which can be added to the other 
even parity Lagrangian pieces. 

In PG, the field strengths T a and R a/3 are untwisted 2- 
forms, similarly the potentials ^ Q and F a/3 are untwisted 
1-forms. One may compare the case of electrodynamics 
with the untwisted potential A and the untwisted field- 
strength F — dA (the differential d is untwisted). Con- 
sequently, a twisted Lagrangian, according to (19), leads 
to the excitations H a and H a p being twisted 2-forms and 
the material currents S a and T a p, see (20), being twisted 
3- forms. 

The Einstein-Cartan Lagrangian 

^ EO = ^^A*(^A^) = ^ii (39) 

£K AK 

is twisted, since the Hodge star in our formalism (see 
[16], Sec. C.2.8) is twisted, that is, it maps twisted into 
untwisted forms and vice versa. The Maxwell Lagrangian 

V Ma , = -^F A*F (40) 

is also twisted, hence of even parity, where Yq is the 
(scalar) vacuum admittance. We recognize that an odd 
number of stars occurring in a Lagrangian, which is ex- 
pressed in terms of field strength and potentials, guaran- 
tees its standard twisted nature. In contrast, the topo- 
logical Chern type Lagrangian 

W = -y^AF (41) 

is only twisted (even parity), if we declare Y\ to be a 
pseudoscalar; the analogous is true in gravity for 

V EC ' = ^jR a0 ^i^A^) = ^-*X, (42) 

AK ZiK 

with the pseudoscalar constant k'; note that this La- 
grangian vanishes identically in a Ricmannian space, 
since = 0. 

D. Parity violating admixtures to the 
Einstein-Cartan Lagrangian 

Already in 1964, Leitner and Okubo [4 !] wondered 
about possible odd parity terms in the gravitational La- 
grangian. Related questions were addressed by Hayashi 
[44] and Hari Dass [45]. The effect of adding the non- 
Riemannian odd parity pseudoscalar curvature to the 
Hilbert-Einstein-Cartan scalar curvature was first stud- 
ied by Hojman, Mukku, and Sayed [46] (for Mukku's re- 
cent view see [47]): 

^hms = ^-(a *R + b *X) 
2k 
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= -^(«o^+W^)A^. (43) 

Note that on the right-hand-side of this equation the star 
only enters in r\ a p, = *(0 a A -dp), that is, the second 
term aj 3 = $ a A tip is of odd parity and thus bo is a 
pseudoscalar. The excitations turn out to be 

H a = 0, H aJ3 = ^-{a O 7] a0 + b tf a (s) . (44) 

ZK 

If we introduce the left-hand-side of (28) as the Einstein 
3- form G a := \r\ a p 1 A i?^ 7 , the field equations read 

a G a — b *X a = kE q , (45) 

yT"A r) a ^ + b T [a A fyj = KT aP . (46) 

The situation with respect to the second field equation 
(46) is similar as in the EC-theory. For vanishing mate- 
rial spin, r Q( 3 = 0, the torsion vanishes, too, T a = 0. This 
can be shown by substituting the irreducible decompo- 
sition of torsion (9) into (46) and using the geometrical 
identities 

(1) T^A Va ^ = 2*«T [a A fyj, 

{3) T»A Va ^ = -^T [a A 00. (47) 

To supply non-vanishing torsion we either need material 
spin or at least (for the vacuum case) field Lagrangians 
quadratic in the field strengths. In this sense, the La- 
grangian (43) is as degenerate as the one of EC and it is 
natural to turn to quadratic odd parity Lagrangians. 

Subsequently questions related to the Vhms - 
Lagrangian, in the realm of classical Riemann-Cartan 
spacetime, were investigated by Nelson [48], Nieh and 
Yan [49] (see also Nieh's recent article [50]), and McCrea 
et al. [51, 52], see also Refs. [53, 54]. 

In the general context of the Ashtekar formalism or, 
more generally, of loop quantum gravity, compare Kiefer 
[55] and Rovelli [56], the VnMS-Lagrangian was taken up 
again by Hoist [57], Freidel et al. [58, 59], Khriplovich et 
al. [60], and Bojowald et al. [61]; similar parity violating 
pieces were studied by Mukhopadhyaya et al. [62, 63], see 
also the related papers by Mielke [64, 65]. 

In the framework of a quantum field theoretical con- 
text, Poplawski [' ], Randono [67], and Bjorken [68] de- 
veloped cosmological models with torsion and parity vi- 
olating pieces that are induced by the vacuum structure. 

Jackiw and Pi [69] proposed a specific model with vio- 
lation of parity and Lorentz invariance in the context of 
GR. They introduced, in addition to the Hilbert-Einstein 
Lagrangian, an external scalar field /3, not to be varied 
in the action principle, multiplied by the Chern-Simons 
(CS) term attached to the curvature: 

VbRCS = ^-*R+ A Rp . (48) 



This model was extended to EC-theory by Cantcheff [ ] , 

Veccs = ^-*R+^-RJ /\Rp a , (49) 
2k 2g 

see also Ertem [ ]. Of course, both theories differ in 
their physical content, as does GR from the EC-theory. 
We turn our attention to the EC-version in (49). 

We know from geometry that a CS-term is an exact 
form. We have for a RC-space, see [15], Eqs.(3.9.3) and 
(3.9.8), 

- ^Rj A R p a = dC KK (50) 

with 

Crr := -^(Vj ARp a + hj AT^ AT 7 a ). (51) 
Thus, 

^eccs = ^~*R- -dCRR. (52) 
2k g 

This Lagrangian contains an odd parity piece quadratic 
in curvature, see (50). However, it is of a fairly degener- 
ate character. Still, since f3 is a prescribed field, the field 
equations are affected by the CS-term. We will come 
back to an explicit evaluation of the curvature square 
piece below. 

Explicit odd parity curvature square pieces were 
present in some cosmological models with spin + and 
spin 0~ modes. These models motivated us for a further 
search in this direction. 



E. Interlude: the impact of the cosmological model 
of Shie-Nester-Yo on PG 

Shie, Nester, and Yo [11, 12], in the framework of PG 
and in accordance with [39, 40], formulated a new cosmo- 
logical model. It contains, besides the graviton mode 2 + 
of GR, one propagating connection and accordingly one 
propagating torsion mode of spin + ; the plus + refers to 
the positive parity. The corresponding Lagrangian reads 
effectively 3 

V SNY = ^ (a *R + ia 2 V A *V ) - ■ (53) 

This model has fairly realistic features and encouraged 
further developments. Li et al. [72, 73], following Ref. 



3 Their actual Lagrangian [ !], Eq. (18) contains also squared 
pieces of the axial torsion ( 3 )T a and the tensor torsion WT a . 
However, they find from the second field equation (for vanishing 
spin) that Mr" = ^T" = 0. Accordingly, in the SNY-model 
only the vector piece ( 2 )T a of the torsion is active explicitly. An 
analogous remark applies to (54). 
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[12], investigated the cosmological evolution of the SNY- 
model with advanced numerical techniques. 

In order to embrace additionally connection modes 
with spin _ , that is with odd parity, in a next step, 
Chen et al.[13] generalized the SNY-Lagrangian to 

y SNY / = - 1 - (a *R - 2A f?) 
In, 

+7T- (02 V A *V-a 3 AA *A) 

OK 

~^(w 6 R 2 ~w 3 X 2 ) V . (54) 

It contains, as odd parity terms, the axial vector torsion 
A and the pseudoscalar curvature X, see (8) and (14), 
respectively. However, these odd parity terms are con- 
cealed in an even parity Lagrangian since each of the 
terms A A *A and X 2 rj = X *X contain one explicit star, 
respectively. Accordingly, this model with propagating 
modes 2 + and is of even parity, but contains con- 
cealed the odd parity terms A and X. 

It is then tempting to provide further add-ons, namely 
the parity violating mixed terms A A *V and X *R — 
R*X. This provides a further extension of the SNY- 
model with respective metric and connection propagating 
modes of 2 + and ± . 



IV. THE SHADOW WORLD OF QUADRATIC 
PG-LAGRANGIANS WITH ODD PARITY 
TERMS 

A. Constructing systematically quadratic odd 
parity Lagrangians 

As discussed by Sozzi [42], from the validity of the 
CPT symmetry in nature and from the fact that the C 
and CP symmetries are only broken in weak interaction, 
but are valid otherwise, one would expect roughly an 
equal amount of matter and antimatter in the Universe. 
It appears however, as shown by Steigman [75] , that only 
what we call matter is around in the Universe. Accord- 
ingly, we have to face the matter/antimatter asymmetry 
in nature and may want to approach this question from 
a gravitational point of view. Can we, in PG, construct 
odd parity Lagrangians in a natural way so that we can 
estimate the possible influence of those terms for the evo- 
lution of the Universe? 4 

We proved identities of the type (32) under heavy use 
of the computer algebra system Reduce, including the 



4 At the Large Hadron Collider (LHC) in Geneva the detector 
LHCb (the 'b' = beauty refers to the 'bottom' quark) was con- 
structed particularly for investigations of matter-antimatter in- 
teractions, namely for the CP-violation in the interaction of b- 
hadrons. It is hoped that these experiments will provide new 
insight into the interrelationship between matter and antimat- 
ter. 



package Excalc for handling directly exterior differen- 
tial forms, see [76-79]. In this way, we can construct 
a quadratic odd parity Lagrangian that in its general 
structure reflects the even parity Lagrangian V+ in (30). 
It is the 'shadow' of V + : 

V_ = -^^R aP A^ 
In 

+1 (a^T a A «T Q + a 2 ^T a A < 3 >T„) 

(»i W R a0 A + ^R^ A ^R a p 

+ ^R a P A ^R a0 + ^R a0 A (b) R a ^ . (55) 

All its constants are pseudoscalars. In a Riemannian 
spacetime such an extended shadow does not exists, 
since only the /^i-term survives, because of T a = and 
(2) Ra p = (3) Ra p = (5) Ra p = q Therefore, in PG we can 

free ourselves from the constraint to use even parity La- 
grangians with only one odd term; in fact, PG brings the 
existence of numerous odd parity Lagrangians to light. 

The special cases of Jackiw & Pi [69] and Cantcheff [70] 
can now be straightforwardly evaluated. If we define in 
four dimensions the (pseudo-) scalar product (A, B) for 
any two tensor-valued two-forms A ai a and B ai a 
by 

(A,B) := *(A ai ... ar A B ai "' ar ) , (56) 

we can write the curvature square piece of our odd La- 
grangian (55) as 

On the other hand, the ECCS-Lagrangian reads 

Vkccs = ^*R+^-*(R,R)- (58) 
2k 2g 

If we substitute the irreducible pieces of the curvature 
into the scalar product, we find [see [15], Eq.(B.4.37)] 

(R,R) = ^R)+2(WR, MR) 

+2(( 3 >i?, ( ^R) + (('%, ^R). (59) 

In a Riemannian space of the theory of Jackiw & Pi only 
the conformally invariant Weyl square piece, the first 
term on the right-hand-side of (59), is left over, whereas 
the RC-space complicates the structures by additional 
post-Riemannian pieces. Comparing (57) with (59), we 
find 

curv 2 I 

V_ (/H = 1, /i 2 = 2, fx 3 = 2, M4 = 1) = —dC RR . (60) 

Q 

In other words, our odd parity curvature square La- 
grangians, for the coupling constants specified, becomes 
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an exact form. Consequently only three of the four /z's 
can be chosen independently. 

The mixed quadratic Lagrangian with even and odd 
parity is then 



v± = v+ + V. , 



(61) 



see also [14]. 



B. Cosmological model with parity violating terms 

If we compare V± with Vg NY ', the next step of 'min- 
imally' generalizing (54) and hopefully keeping the nice 
properties of the model is to allow for unconcealed odd 
parity pieces, but only those odd parity pieces that al- 
ready occur in (54), namely X and A. Thus, starting 
with (CI) and putting the following constants to zero, 

a 1 = 0, <j\ = 0; Wi — w 2 = W4 = w 5 = 0; 

Mi = A*2 = fH = 0, (62) 

we arrive at the new Lagrangian 

Kbhn = i (-a (6) i? Q ' 3 A Va p - b {3) R afl A tf Q(3 - 2\ n 
+a 2 {2) T a A * (2) T a + a 3 (3) T a A * (3) T Q 
+2a 2 (2) T a A (3) T Q ) 

i (w 3 ^R a ? A *WR af) + w 6 ^R afi A <^R af) 



+^R afi A^R a p 



(63) 



Substituting the irreducible pieces of the torsion (7) and 
(8) and of the curvature (14) and (15) into (63), we 
find the more compact form of our new gravitational La- 
grangian 

Vbhn = ~(a *R + b *X - 2A t?) 

ZK 

+— (a 2 V A *V - a 3 A A *A - 2a 2 V A A) 
6k 



24g 



(w 6 R *R - w 3 X *X + fi 3 R*X) . (64) 



The constants 60, o~ 2 , M3 are psetirfoscalar, the remaining 
ones are scalar. 5 

We can read off from the Lagrangian a symmetry be- 
tween R and X, between wq and 103, and between ao and 
60 ■ There is also a symmetry between V and A on the 
one side and R and X on the other side; this implies that 
02, a 3 j 02 are mirrored in wq,w 3 ,[jl 3 . These symmetries 
are also reflected in the field equations. We will come 
back to this in Sec. IV E. 



A teleparallel Lagrangian with a term with <72 7^ has been 
considered earlier by Miiller-Hoissen and Nitsch [80]. 



The decomposition of the linear terms in R and X into 
Riemannian and post-Riemannian pieces, modulo surface 
terms, yields 



Vbhn 



1 r 

2k 
2 



«o 



*R - 2X n + «T Q A (6 + a T a 



m+VA *V+ -m~AA A - 2 m x V A *A 
3 6 3 



1 

24^ 



{w 6 R *R-w 3 X*X + f i 3 R *X) 



(65) 



with 



a 2 

a a - y i m 



a - 2a 3 , m x := b Q + a 2 . (66) 



We shall see that the m's play a role in the discussion of 
the second field equation. Note that m + and m~ are of 
even parity, whereas m x is odd. A corresponding decom- 
position of the quadratic curvature terms in (65) doesn't 
seem to provide new insight. 



C. Diagonalization of the BHN-Lagrangian 

1. Eigenvalues of the kinetic matrix of the translational 
gauge potential i9 Q 

In the Lagrangian (65) the term V A *A = A A *V = 
A^V^n represents an interaction term of two four- vectors 
of the type vector x axial vector. To get some more 
insight into the dynamics of the fields governed by the 
Lagrangian (65), we decompose these four- vectors into 
(18 3) with 



Aa=(A ,A) and V„= Vo.V 



(67) 



The translational part of the Lagrangian (65), with 
(l) rpa _ j g proportional to the quadratic form 

Q := 4(V 2 - V 2 - Vf - V|)m+ 
-(Al-Al-A 2 2 ~A 2 3 )m- 
+4{A V - AM - A 2 V 2 - A 3 V 3 )m x . (68) 

Expressed in terms of matrices, we have 



(V ,V,A ,A) 



— 4m + g — 2m x g 
— 2m x g m~g 



(V ,V ,A Q ,A) - T ■ (V ,V ,A Q ,A) T , 



V 
A 
\AJ 



with g as the four-dimensional Minkowski metric. 
With the useful abbreviations 



x : = 



4to + . 



V 



2m y 



m 



(69) 



(70) 



10 



the new matrix T reads 

fx 

-x 




V o 

-y 


\o 

xg -yg 

yg zg 



T 



o 



x 



-x 

-z 





£/ 






v o 





-v o 





z 





o \ 


// o 
y 



(71) 



It has some simple properties: 

r T = r, 

trace T = — 2(x — z) , 



detT = 

T- 1 = 



(xz + y 2 ) 



1 



>0, 

~zg -yg 
-yg xg 



(72) 



In the following we will make use of Dirac's bra-ket 
notation, see Dirac [81] and Schouten [82]. For this pur- 
poses we define (abstract) vectors according to 

<X\ := (V ,V,A ,1) and \X >:= (V , V , A , A ) T 

(73) 

such that the quadratic form Q becomes 



Q =< X\T\X > 



(74) 



To diagonalize the form (74), we introduce a new vector 
\y > together with a suitable orthonormal matrix K, such 
that 



\X>=:K\y> and <X\=<y\JC I 



(75) 



Substitution of (75) into (74) yields the covariant expres- 
sion 

Q = < X\T\X >=<y\K T -r -K\y > 
= < y\v\y > . (76) 

We will choose the orthonormal matrix JC such that the 
product T> := JC T -T-JC is diagonal and thus contains the 
eigenvalues of T as entries. 

The eigenvalues of T turn out to be: 

Ai = - (x - z + ^J(x + z) 2 + 4y 2 ) , 
A 2 = A 3 = A 4 = -Ai , 



A, = 



- (x - z - ^{x + z) 2 + 4y 2 ) , 



A 6 = A 7 = A 8 = -A 5 



(77) 



According to the Lorentz structure, that is, the (1 © 3)- 
decomposition of four-vectors, we have two different 
eigenvalues for the time components and 2 threefold 
eigenvalues for the spatial components, respectively. 



For the explicit construction of the matrix K, we need 
the eigenvectors u\ n of T, where n runs from 1 to 8. 
Their components can be expressed (for y ^ 0) in terms 
of A and B : 



A: 



— (x + z + y/(x + z) 2 + 4y 2 ) , 

B:=-— (x + z - y/(x + z) 2 + Ay 2 ) . (78) 
2y V I 

We normalize the eigenvectors by means of 

a := y/l + A 2 and b := \J\ + B 2 . (79) 

The columns of the matrix K, are the normalized eigen- 
vectors of the matrix T. Accordingly, we find K. = 



(- 


-A/a 











-B/b 








\ 





1/6 











-l/a 
















-A/a 











-B/b 
















1/6 











-l/a 




If a 











1/6 
















B/b 











-A/a 
















l/a 











1/6 
















B/b 











-A/a) 



(80) 

Because of the identity BA = — 1 and its consequences 
—A/a = 1/6 and B/b — l/a, there is actually much more 
symmetry than is apparent in the matrix (80), which in 
fact depends essentially on only one parameter. This 
matrix has simple properties like det/C = +1, JC ■ /C T = 
1-8x8; the eigenvalues are e , both with multiplicity 
four and with tana = 1/A, cf. Eq. (87). Hence in the 
eight-dimensional vector space the matrix K, represents a 
pure rotational matrix. 

Thus, the quadratic form Q assumes the diagonal form 



Q =< X\T\X > 
=< y|diag(A x ,- 



A U A 5 ,-A 5 )\y > : 



(81) 



with the obvious abbreviations Ai := A], (1,1, 1) and 
A5 := A5(l,l,l). The new vector \y > turns out to 
be 



|y>:= 



V 

A 



= JC T \X >=fC T 



V 
A 
\AJ 



(82) 



Accordingly, the quadratic form (68) can now be written 
in diagonal form as 



Q = Ai ( V - V 1 



A, An" - A" 



(83) 



We need to pay special attention also to the case y = 0, 
namely when there is no coupling between the vector and 
the axial vector of the torsion. In this case, the matrix 
T is diagonal with 



T(y = 0) = 



-xg 
zs 



(84) 
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and the eigenvectors of (84) correspond to the eight- 
dimensional unit vector. Thus, for y — > 0, the trans- 
formation matrix becomes the unit matrix: K, = lsx8- 
The transformation to the diagonal matrix yields 



/C T T/C = T> = diag(x, —x, —x, —x, —z, z, z, z). 
In this case, the quadratic form Q reduces to 



Q(y^0) = :r(V 2 -V 2 ) - z ( A 2 - A 2 



4m+ V -V 2 -m" ( A - A 



(85) 



(86) 



2. Representation of K, in terms of angles 

The matrix K. as a rotation matrix can be parametrized 
by introducing a suitable angle variable. For this purpose 
we use 



1 



sin a 



VI + A 2 



and cos a 



A 



(87) 



In terms of this parameter a, the matrix JC can be dis- 
played in compact form as 



cosal4x4 sinal4x4 
— sinal4 X 4 cosal4 X 4 



where I4X4 is the 4x4 unit matrix and (x, y, z) are related 
to the parameter a by 



tana 



x + z — j3 



-2y 



X + z + p ' 



with 



13 := v/(x + z) 2 + iy 2 . 
Thus, the matrix product T ■ JC reduces to 



(89) 



(90) 



T-JC=( cosaI ^ shiaI ***).V = lC-V, (91) 
y — smal4x4 cosal4x4 J 

the condition for a similarity transformation. 

3. A possible parameter set 

We can read off from (77) the relations A1A2 < and 
A 5 A 6 < 0. As an example, let us consider the case A x > 
and A5 > 0. With these assumptions we can define 

i/:=\/AiV, Ai > 0, 
a*:=VA 5 A, A 5 > 0. (92) 



^ :=vAiV o 
a :=VA 5 A 



With (92), the quadratic form Q (83) assumes its 'special 
relativistic' appearance 



Q 



v 2 ) 



an 



a 



-9^ {v»v v + ^u v ) 



(93) 



This is the sum of the squares of two four-vectors in a 
suitable orthonormal reference frame and the Lorentz co- 
variance is manifest. 

Let us analyze the conditions to be imposed provided 
one assumes Ai > 0, A 5 > 0. From (77) we derive the 
constraints 



z > 



Am — m > 

3a - 2(a 2 - a 3 ) > , (94) 



and 



xz + y 2 <0 ^==^ m + m +(m x ) 2 <0 
(ao - y ) (oq - 2a 3 ) + (b + a 2 ) 2 < . 



(95) 



On the other hand, assuming instead Ai > 0, A5 < 
leads to the condition 



y 2 >o 



m m + (to x ) > 



[a - y ) (ao - 2a 3 ) + (&o + cr 2 ) 2 > . 



(96) 



One could similarly find the parameter conditions asso- 
ciated with the other two cases. 

If and A p arc both timelike — as they will turn out 
to be for the cosmological model which we derive below — 
every set of parameters fulfilling the inequalities (94) and 
(95) will lead to a strictly positive kinetic energy matrix 
for the translational gauge potentials. 



4- Eigenvalues of the kinetic matrix of the Lorentz gauge 
potential Y al3 



Similar as in Sec. IV CI, we consider the quadratic 
form C representing the curvature square terms in (65) 
and hence the kinetic parts of the connection. This 
quadratic form is given by 



C 



:= w 6 R 2 - 
= (R,X) 
= <Z\B\Z> 



w 3 X 2 + fj, 3 RX 



( w 6 


AW2\ , 




Us/2 


-w 3 J 





(97) 



We wish to diagonalize the symmetric (2 x 2)-matrix 



B 



w 6 fj, 3 /2 
H 3 /2 -w 3 



(98) 



whose eigenvalues are 
Ai,2 = ~w 3 ± y/x), x •= {w6 ~ W3) 2 + A, (99) 



where 



A 



-4deti3 = 4w 3 wq + /i 3 . 



(100) 



This diagonalization can be simply accomplished with a 
rotation matrix whose columns are orthogonal unit eigen- 
vectors; the transformation matrix and the eigenvectors 
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can be parameterized by a single angle 7 that can be 
determined from 



tan 7 = 
Then, 
M T BM = 



w 6 + w 3 



'X 



-H3 



/'3 



Ai 
A 2 



w 6 + w 3 + y/x 



M := 



cos 7 sin 7 
— sin 7 cos 7 



(101) 



(102) 



Consequently, with 

(R,X) = (R,X) ■ M, (103) 

the quadratic form C can be expressed as 

C = < Z\B\Z >=< ZM\M T BM\M T Z > 
= <Z\D R \Z>= X 1 R 2 + X 2 X 2 , (104) 

where T>r := A4 T BA4 denotes a diagonal matrix. 

The quadratic curvature terms in the Lagrangian can 
now be rewritten in the form 



Vt 



R 2 



24g 



R 2 



X2X 2 



(105) 



As an illustrative example, consider in particular the 
case where Ai and A2 are both negative. This immedi- 
ately leads to the following constraints on the coupling 
constants: 

^6 — "^3 < , and /ig + Aw 3 w§ < , (106) 

from which we can infer that 

w e < 0, and w 3 > (107) 

for the aforementioned case, Ai < and A2 < 0. Then 
we can rescale the variables and introduce new ones ac- 
cording to 



R 



R 



and X := 



X 



(108) 



such that the quadratic form C in this particular case 
becomes 



C = - (R 2 + X 2 ) 



(109) 



For the three other cases, namely (Ai > 0, A2 > 0), (Ai < 
0, A 2 > 0), (Ax > 0, A 2 < 0), one can do an analogous 
rescaling. 



5. Partly diagonalized Lagrangian 

Thus, the process of diagonalization for the case Ai > 
0, A 5 > 0, Ai < 0, A2 < leads to the following diagonal 
pieces of the VeHN-Lagrangian, 



^[M-fVW-a-H-i.jLa,, 



Vr 2 = WT- R 
R 2Ag v 



(110) 



(with analogous results for the other sign choice cases). 
Collecting the results received so far, we can give a new 
representation of the Lagrangian (65) in the form of 



Vt 



BHN 



= i-(a 5-2Ao) ?7 

+ fao * {1) T a A {1) T a + b {1) T a A w T a ) 

+ ^Qri-^ r C? 1 . (Ill) 
12k 24g 



6. Correspondences of eigenvalues of the kinetic matrices 
to spin & parity 

In this section we considered sufficient conditions for 
the coefficients of the kinetic energy matrix being posi- 
tive. We now assume that the trace and the axial-vector 
pieces of the torsion are both propagating independently. 
Then we require that the four- vectors and a M are time- 
like during the whole evolution. Accordingly, the propa- 
gation of independent massive modes is characterized by 
v^v^ < and a^a^ < 0. This is met by the requirement 
(92). Other choices of the signs of A^'s will lead to space- 
like four-vectors. The null case will be treated separately 
in a continuation of this paper. 

The Lagrangian (111) admits the introduction of a 
number of strictly positive functions, that is, functions 
of Afc 's, such that for each of its dynamical variables we 
can associate to each eigenvalue of the kinetic energy ma- 
trix the corresponding spin & parity state. 

If we decompose the four-dimensional one-forms V and 
A into (1 © 3), respectively, we are naturally led to their 
spin contents. Namely, we can introduce (massive) three- 
dimensional vectors v and a and can associate to each of 
them a corresponding three-dimensional spin & parity 
state. This method is not sensitive to a possible occur- 
rence of multiplicities of spin & parity states. For this 
purposes, we have to investigate the corresponding lower- 
dimensional subspaces. 

In our model, with ^T a — 0, we have only propa- 
gating scalar and three-dimensional vector modes. In a 
heuristic manner, the diagonalization allows for the fol- 
lowing tentative correspondences: 



Ai(Vo) 

Ai(V) 
A 5 (A ) 

A 5 (A) 
Xi(R) 
A 2 (A) 
R 



+ 
1+ 
0" 

1" 

0+ 
0" 
2+ 



(112) 



In the case of a non- vanishing tensor torsion ^T a 7^ 0, 
those terms deliver massive modes of spin state 2 that 
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would combine with the corresponding spin 2 + -mode of 
the Riemannian curvature scalar R. Because of the com- 
plexity of these results, we will defer their presentation 
to follow up work. 



D. Excitations of the gravitational held 

We differentiate the Lagrangian (64) with respect to 
torsion and curvature. Then, with the help of (19), we 
find the translational excitation 

-^ 2 ( (2) T Q +( 3 )T Q ) 



(113) 



and the Lorentz excitation 
1 



H. 



a/3 



o v 



2 k 
1 

Q 



+Y e ^ ( {3) R a p + (6) R a p) , (H4) 

respectively. Alternatively, the field excitations can be 
given in terms of the field strengths in a more suitable 



and symmetric form (V = V a $ a and A = A a $ a ) 

H a = — [{a 2 V l _ l - (T 2 A tl ) rf a 

+ (^ + 03^)^], (H5) 



From (115) we find in particular 



H a A zT = ~*(a 2 V-a 2 A) 

K 



(117) 



E. Explicit form of field equations for the new 
Lagrangian 

By substituting the excitations (113), (114), and the 
Lagrangian (64) into the gauge currents (23) and (24), 
and then the latter two, together with the excitations 
(113), (114), into the field equations (21) and (22), we 
find the explicit forms of the field equations: The first 
field equation reads, 



T^R - T7T X G « H Va 

k bg lzg J k 



k bg lZg 



+ ^(v.v-iv«)* + ^(^-ix'«)* 

+i (<=«-. W T") A [«, + m «< 3 >J>+ n (< 2 >7i + <»>3>)] = £„ , 
and the second field equation, 

^2 * (1) T [Q A 0fl - I V A tic/, + 1-AA Kp) + ^ (2 (1) T [Q A 0^ - |v A d a p Va t 



2 k 

+ 2Jg ( 2w3dX ~ ^ dR ) A ^ - 24" {2w&dR + ^ dX ^ A ~ 24" ( 2wqR + T1 A 

+ ]Jg ( 2w3X ~ /i3i? ) T I Q A ' d/3] ~ Jk [ a2V [«^] ~~ °2.A[ Q J7/j] + (o 3 -4 + cr 2 V) A Q(3 ] = r a , 3 . 



(118) 



(119) 



r 



The source of the second field equation is the material 
spin angular momentum 3- form r a p. According to its 
definition in (20), it is antisymmetric in a and /3. It 
is related to the source of the first field equation, the 
canonical energy-momentum 3- form of matter E a , via 
the angular momentum law DT a p + 0r a A E^i = 0. 



It may not be superfluous to look at the struc- 
tures of the two field equations and at the ways we 
ordered them. The first line of (118) emerges from 
the curvature dependent pieces of the gauge energy- 
momentum E a . We find, symbolically written, ~ 
Einstein + cosmol. term + pseudoscalar curv + curv 2 . 
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If only ao ^ and bo ^ 0, we recover the left-hand-side 
of (45); if only a ^ 0, we have just the EC-theory. 
The second and third lines of (118) are of the following 
structure: d torsion + torsion 2 . From the point of view 
of gauge theory, 'd torsion' is the leading term, see DH a 
in (21); the remaining 'torsion 2 ' pieces collect the torsion 
dependent parts of the gauge energy- momentum E a . Of 
course, also the frame e Q and the coframe *& a feature in 
this equation directly or indirectly via rj a = *d a . 

We ordered the second field equation (119) in a similar 
way. In the first line we have the terms linear in torsion 
originating from the gauge spin E a p, see (24); compare 
also as special case the left-hand-side of (46). In the 
second and third line, we have 'd curv + curv x torsion'. 
The leading term is 'd curv', the rest arises from the 
differentiation process of 'D curv' with the help of (2). 
Again the coframe $ a enters explicitly or implicitly via 
= $a A &p, r] al3 = *(tf a A dp), and 77 a/ 3 7 = *("d a A 
A0 7 ). 

Accordingly, the two field equations are now expressed 
in terms of the torsion and the curvature of spacetime. 
In the sense of gauge theory one may now want to in- 
sert the definitions of torsion, Eq. (5), and curvature, 
Eq. (12) cum (11). Then we would get second or- 
der quasilinear partial differential equations (PDEs) in 
coframe and Lorentz-connection: i d*dd+ lower order ~ 
energy-momentum' and l d*dT + lower order ~ spin'. 
However, our experience on the search for exact solu- 
tions, in particular for cosmological models, has shown 
that it is to be preferred to stay with the well-behaved 
tensor-valued two-forms of torsion and curvature and 
not to switch to the proper gauge variables coframe and 
Lorentz-connection. 

As we saw already, the Lagrangian (64) has a remark- 
able symmetry which we will find also on the level of 
the excitations (115) and (116), the field equations (118) 
and (119), as well as on the level of the coupling con- 
stants. Hence we can introduce the following tentative 
correspondences (in four dimensions) between variables 
and parameters, respectively, 



Vo- 



A 



«3 



(120) 



Ro— • X 




Wq o • w 3 


G«o— • *X a 




a o — • b a 



(121) 



2 °- 



M3 



(122) 



F. A consequence of a topological term 

There is a subtlety present in the Lagrangian (64) 
and the corresponding field equations (118) and (119). 
Because of the Nieh-Yan identity [49] — see also [15], 
Eqs.(3.9.7) and (B.4.15)— we have 

d($ a AT a ) = T a AT a + R a/3 A$ afl = T a A T a — *X . (123) 

The torsion square term can be expressed in its irre- 
ducible components according to 

T a A T a = (1) T Q A {1) T a +2 (2) T Q A (3) T Q 

= WT a A (1) T a - \aa*V. (124) 

We substitute (124) into the right-hand-side of the Nieh- 
Yan identity (123) and find 

d(&* A T a ) = {1) T a A (1) T Q - \A A *V - *X . (125) 

o 

For the sake of a neater argument let us first extend 
our parity mixed PG Lagrangian (64) by including also 
the CTi-term from (55): 

V = 1/bhn + V= Vbhn +— {l) T a A V>T a . (126) 



If we substitute (64) into (126), then we recognize that we 
can recover from kV the right-hand-side of the Nieh-Yan 
identity (125) for the specific coupling constants 



b = -2, o\ = \ 



(T 2 =2, 



(127) 



all other constants, apart from k and g, vanish. Since 
the left-hand-side of (125) is an exact form, the choice 
(127) corresponds to a 'null Lagrangian' with vanishing 
field equations. 

By the same token, we can add a multiple (say e/fc) of 
the exact 'topological' form d(d a A T a ) to V. After some 
simple algebra, we find 



V + -d{$ a A T a ) = -(tr x + e) (1) T a A (1) T Q 

K 

[ao*R + (b Q - 2eyX - 2Xo V ] 



1 

2k 
1 

6 k 
1 

24g 



[a 2 V A *V — a 3 A A* A- 2(a 2 + 2e)A A *V] 
(w 6 R*R-w 3 X*X + fi 3 R*X) . (128) 



This is equivalent to certain changes in the parameters 
of our action (126), specifically 

oi ax + e, b Q -> bo - 2e, a 2 ->• a 2 + 2e. (129) 

From this we can infer that the field equations can- 
not depend on the parameters a\, bo, a 2 by them- 
selves, but rather must depend on these parameters only 
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through certain combinations which are invariant under 
this transformation, such as 



b + 2(7 1 , 2cti - er 2 , 



(130) 



The two field equations of the Lagrangian V are found 



via H a = —1o\ ^Ta/n as 



l.h.s. of (118) 



2(7l 



D {1) T a 



l.h.s. of (119) - — % A Wt 0] = r Q/ 
Specifically, using the first Bianchi identity, 



(131) 



(132) 



*X a = R 0a A^e DT a = D(«T Q + &T a + ^T a ) , 

(133) 

and the expressions (7) and (8) for the irreducible pieces, 
one can find that the bo, a±, and 02 terms on the left- 
hand-side of the first field equation (131) add up to 

~D [(b + 2a 1 )«T Q + (b + o 2 ) ^T a + ^T a J 

+i(2a t - a 2 )(e aJ «T^) A (V>T P + ( 3 >7>) ; (134) 

similarly, for the left-hand-side of the second field equa- 
tion (132) we have 

- (6 + fcri) (1) T [Q A 0fl - — (6q + o- 2 ) (2V A i) a0 +A A r? Q/3 ) . 

(135) 

There are several points worth noting: (i) As expected, 
the parameters occur only in the invariant combinations 
(130). (ii) All the o\ terms are proportional to ^T^. 
(iii) For solutions such that Wt^ 1 vanishes, the equations 
contain the parameters only in the combination &o + <T 2- 

Having obtained these insights, it is no longer neces- 
sary to keep the rather complicated <j\ term. Exploit- 
ing the freedom to choose a suitable e in (129), namely 
e = — (7i, we can, without loss of generality, drop the a\- 
term altogether and return to our model Lagrangian (64) 
with the two field equations (118) and (119). 



V. FRIEDMAN COSMOLOGIES WITH 
PROPAGATING MODES OF THE 
LORENTZ-CONNECTION 

Since the early 1970s cosmological models for EC and 
PG have been developed, see Kopczyhski [87], Trautman 
[88], Tafel [89], Kuchowicz [90, 91], Kerlick [92], and oth- 
ers [93, 94], to name a few. Minkevich et al. [95] devel- 
oped the subject in a series of papers. A report on the 
status of the subject was given by Puetzfeld [96]. 

For our new Lagrangian we follow these procedures 
and search for FLRW type of cosmological models. 



A. Homogeneous and isotropic coframe and torsion 

Assuming a homogeneous and isotropic scenario, the 
orthonormal coframe for a Friedman cosmos is 



if 

t? 2 
tf 3 



dt 

a(t)dr 
7! - kr 2 
a(t)rd9 
a(t)rsin 6d<j> 



(136) 



with the metric 



a=l 



-dt 2 



o 2 (t) 



dr 2 +a 2 (t)r 2 (d9 2 + 



sm 



(137) 



where a(t) is the expansion factor and k the curvature 
index. 

Now we can compute, up to antisymmetry, the non- 
vanishing components of the Ricmannian connection 
(a, b,c, ... = 1,2,3 are spatial anholonomic (frame) in- 
dices): 



rV = 

fa 1 = 
f 3 2 = 



a(t) 

a(t)r 
VI - kr 2 q3 

T\ $ 

a(t)r 
cot 9 „o 



a(t)r 



(138) 



Like in Einstein's theory, the temporal rate of change of 
the expansion factor a(t) determines the Hubble function 

H(t) := a'(t)/a(t) . (139) 

By differentiation and elimination of a'(t), we find 

H'(t) + H 2 (t) = a"(t)/a(t). (140) 

This determines the Riemannian sector of spacetime. 

The most general torsion compatible with homogeneity 
and isotropy can be characterized by two independent 
functions u(t) and v(t), see Goenner & Muller-Hoissen 
[83] and Baekler [84]. We will choose for the torsion the 
parametrization 



rpQ 

T 1 
T 2 

rp3 





u{t)$ 01 + v(t)$ 23 
u(t)$ 02 + v(t)V 31 
u(t)# 03 + v(t)$ 12 



(141) 
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The irreducible decomposition of T a implies a vanishing 
tensor piece, whereas the vector and axial vector pieces 
survive: 



(1) rpa _ q (2)jia _ 




(3) T « = v 




(142) 



This yields for the corresponding 1-forms in (7) and (8) 

V = -3u(t)tf , A = -3v(t)$° . (143) 

Incidentally, the purely spatial part of the torsion ^T a 
corresponds to Cartan's spiral staircase [85, 86] with a 
time dependent pitch of the spiral. As such, one can 
easily visualize it, see [86]. 

By simple algebra, we can calculate the contortion K a p 
in terms of the torsion, see [15], 



K a/ 3 = e[ a jTp] - - (e a AepjT 1 )'d 1 . 



(144) 



We find then the Riemann-Cartan connection according 
to r Q(9 = T a/ 3 - K a p or, explicitly, 

r a ° = [H(t) - u(t)} $ a | 

T a » = f a b +\v(t)e ca b r ) ■ (M5) 



B. Irreducible pieces of the curvature 

Having now coframe, connection, and torsion at our 
disposal, we can immediately calculate the different 
pieces of the curvature 2-form. We introduce the Hubble 
function H{t) (139) and find straightforwardly 



(!) R 
R a p 

a /3 



± V, 



0, 

hv(t)[H(t)-2u(t)\ 



At)} 





#23 


_tfl3 


012 








_^03 


$02 





o 







U 


o 


o 






i3) R af 3 



- [H'(t) - u'(t) + H(t)u{t) - u 2 {i) 



/0 ti 01 V 02 tf 03 N 

o $ 12 tf 13 

o o d 23 

\o o o 



{5) R a (3 = 0. 



{&) R aP - 

with the (pseudo)-scalar functions 

X(t) = -3{v'(t) + v(t)[m(t)-2u(t)]} , 



(146) 

(147) 
(148) 

i 

(149) 

(150) 
(151) 



(152) 

R{t) = 6{{H'(t)-u'(t)]+H(t)[H(t)-u(t)] 

+ [H(t)-u(t)} 2 - 1 -v 2 (t) + -^ ) Y(153) 

The matrices in (147) and (149) are antisymmetric, re- 
spectively. Therefore, we indicated those matrix elements 
with a diamond symbol that are, because of this antisym- 
metry, redundant. 

Since we chose as our variables curvature and torsion, 
the relations between curvature and torsion provided by 
the nonvanishing irreducible pieces and, in particular, by 
(152) and (153), are vital information in our search for 
exact solutions. 



C. A spinless perfect fluid as model of matter 

The sources of the two field equations are the energy- 
momentum current S Q and the spin current r a p. These 
three-forms we represent as tensor-valued zero-forms ac- 
cording to 

S Q = Ta^Vfi , T af3 = Sap^T]^ ] (154) 

the reciprocal relations read 

rj = * (e q a , s a0 -y = * (r a0 a n) . (155) 

In the following we will only consider matter models 
with vanishing spin, r Q/ g = 0. This simplifying assump- 
tion, which may be justified for spherical symmetry, cer- 
tainly has to be dropped in a more advanced stage of our 
model building. 

Due to the Friedman (or FLRW) symmetry of our 
cosmological model, the energy-momentum tensor must 
have the spinless perfect fluid form 

Tj = [p(t) +p(t)]U a UP+p(t)6P, (156) 

where p = p(t) is the energy density, p = p(t) the pres- 
sure, and U a the four-velocity of the fluid, with the 
normalization U a U a = — 1. Because of the symme- 
try requirements, we only have the dependencies p(t) 
and p(t). In a comoving reference system, we have 



,J o > ' 



(157) 



Tj = diag (~ P (t), P (t),p{t),p(t)) . 
From (156) we deduce for the trace 

T a a = -p(t) + 3p(t) (158) 
and for the trace-free part 

Tt fi = W)+p{t)](u a U' 3 + \5^ 

= 7b(t)+p(t)]diag(-3,l,l,l). (159) 
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D. Differential equations for torsion and curvature 



1. First field equation 



According to the FLRW-symmetry requirements, the 
first field equation (118) as a vector-valued 3-form, has 
only two (algebraically) independent components and in 
the same manner the second field equation (119) as an 
antisymmetric tensor-valued 3-form also has only two in- 
dependent components. 



We substitute into the first field equation (118) the 
coframe (136), the torsion (141), and the expansion fac- 
tor via (139) and (140), but leave R(t) and X(t) as they 
are. Then we find as independent non-vanishing equa- 
tions only the components p(t) — e% j[e2_i(e3_i£o)] & n d 
p(t) = eoj[e 2 j(e 3 jS 1 )]: 



K P (t) = l - [a R(t) + b X(t)} + 3a {[u(t) - H(t)}' + H(t)[u(t) - H(t)}} + h [v'(t) + H(t)v(t)] 



3 
2 

K 

'4g 

K 



K 

~2Aq 



{a 2 u{t) [u(t) - 2H{t)] - a 3 v 2 (t)} + 3cr 2 v(t) [u(t) - H(t)} 
[2w 6 R(t) + nnX(t)] {[«(<) - H(t)}' + H(t)[u(t) - H(t)}} 
\p 3 R(t) - 2w 3 X(t)} [v'(t) + H(t)v{t)] 
[w e R 2 (t) + fi 3 R(t)X(t) - w 3 X 2 (t)] - A 



(160) 



and 



Kp(t) 



--[a R(t)+b X(t)}+a Q 



H'(t)+3H 2 (t) 



2k 

am 



{a + a 2 )u'(t) + ~{2a 2 -b )v'(t) 



-(5a + 2a 2 )H(t)u{t) - ^(5b - 4a 2 )H(t)v(t) + i(4a + a 2 )u 2 {t) + (2b Q - a 2 )u(t)v(t) 



1 



(ao + as)v 2 {t) + A 



12g 

K 

2Ag 



[2w 6 R(t) + n 3 X(t)] 



9k 1 
H'(t) + 3H 2 {t) + - u'(t) - bH(t)u(t) + 2u 2 {t) - -v 2 (t) 

a z (t) I 



[H 3 R(t) - 2w 3 X(t)} [v'(t) + 5H(t)v(t) - 4u(t)v(t)] 



+ 24^ [ WeR2 (<) + ^ R ^ X ^ ~ W3X2 W] 



(161) 



A further relation between the fluid density p(t) and the pressure p(t) can be gained by taking the trace d a A E Q of 
the first field equation or, equivalently, by computing p(t) — 3p(t) from (160) and (161). However, in order to find 
a compact expression, we resolve (152) with respect to v'(t) and (153) with respect to H'(t). If we substitute these 
expressions, we find 6 

K [p{t) - 3p(t)] = a R(t) + (b + <T 2 )X(t) + 3a 2 [u'(t) - u 2 (t) + 3H(t)u(t)] + 3a 3 v 2 (t) - 4A 
= ~(2a -a 2 )R[t) + (b + a 2 )X(t) 

+ ^(4a 3 - a 2 )v 2 (t) + 3a 2 ^ff'(t) + 2H 2 (t) + -J^J - 4A . (162) 



Let us now get back to the first field equation (160) with (161). One strategy is to eliminate the time derivative 
H'(t) of the Hubble function by a suitable linear combination of (160) and (161). Accordingly, we put the linear 



6 That on the right-hand-side only a derivative with respect to the 
torsion can appear can be seen from the general structure of this 
expression: $ a A S a = (p - 3p)n = -d(tf Q A H a ) — 4V — T a A 
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combination {ei_i[e2-i(e3_i£o)] + 3eo-i[e2_i(e3_i£i)] — [p{t) + 3p(i)]} to zero and, isolating the derivatives of the torsion 
functions, we find 

- 3 [a 2 u'(t) - a 2 v'{t)\ = k [p(t) + 3p(t)] + 3H{t) [a 2 u(t) - a 2 v(t)} + a R(t) + b Q X(t) 



+6b v{t) [H(t) - u(i)] + — [2w 3 X(t) - p 3 R(t)} v(t) [H(t) - u(t)} 

2p 

_6o |[^(t)-ti(t)] a -i« a (t) + ^ 
+^ [p 3 X{t) + 2w e R(t)} {[#(*) - u(t)f - \v{t) 2 + 
-2A - ^- [u; 6 i? 2 (t) + Li 3 R(t)X(t) - w 3 X 2 (t)] . 



(163) 



Alternatively, the last relation (163) can also be expressed as 
3 d 



a(t) dt 



{a(t) [a 2 u(t) - a 2 v(t)}} = n [p{t) + 3p(t)] + a R(t) + b X(t) + 6b v{t) [H(t) - u(t)} 

— (2w 3 X(t) - p 3 R(t)) v(t) [H(t) - u(t)} - 6a l[H(t) - u(t)} 2 - \v 2 {t) + 

— [p 3 X(t) + 2w 6 R(t)} {[#(*) - u(t)} 2 - \v{t) 2 + ^^- -^u 
~ [w 6 R 2 (t) + p 3 R{t)X{t) - w 3 X 2 (t)\ . 



(164) 



r 



This equation suggests to impose the interrelationship 

CT2 



u(t) = O v(t) 



with 



rt 2 



(165) 



{v(t)(oT u(t)),R(t),X(t),H(t),p(t),p(t)}. We will defer 
the study of these two alternatives to future work. 



We can also manipulate the first field equation (160) 
between the two torsion pieces as a simple special case; with (161) in a different way in order to arrive at al- 



a further possible choice could be 



u(t) — av(t) = — -r- , c = constant ^ . 
a(t) 



(166) 



gebraic relations. We can resolve (153) and (152) with 
respect to u'(t) and v'(t) and substitute these expressions 
into (160) and (161), respectively. 



Then the left-hand-sides of (163) and (164) vanish and 



After eliminating the derivatives of the torsion we ar- 



we find a purely algebraic equation in the variables rive at 



Kp{t) = -3m + u(t)[2H(t) - u(t)]- -m~v 2 (t) + 3a H 2 {t) + 



-3m x [H(t) - u(t)} v(t) - ^ {p 3 X(t) + 2w 6 R(t)] W) u(t)} 2 \v 2 {t) + 



K 

~2Aq 



\p 3 R(t) - 2w 3 X(t)} v(t)[H(t) - u(t)} 



[w 6 R 2 (t) - w 3 X 2 {t) + p 3 R(t)X(t)] - A 



(167) 



and 



K p(t) = --m + R(t)~-m x X(t)+2(m + 
3 3 

+ hn~V 2 {t) + a j [H(t) - u(t)Y 



a )H'(t) + m x v(t) [u{t) - H(t)) 
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+2(m+ - a ) |ff(f) [2H(t) - u(t)] + - U 2 (i) - ^(i) + 
WW - 2™ 3 X(t)] w(t) [#(*) - u(t)] 
|/i 3 Jr(t) + 2w 6 R(t)} \[H(t) - u{t)] 2 - \v 2 {t) + 



I2g ^ w ^ WJ 4 w a 2 (i) 



/v 



Ki? 2 (t) - ™ 3 X 2 (i) + /i 3 i?(<)X(t)] + A . (168) 

Inspecting the equations there are the following dependencies, 

p(t) = p[a{t),H{t),u{t),v(t),R(t),X(t)} , 

p(t) = p[a(t),H(t),H'(t)Mt)At)>R(t)>X(t)] > (169) 

that is, the density depends only algebraically on the variables whereas the pressure, besides non linear terms, contains 
only the derivative of H(t). Thus, also this general case belongs to the class of descriptor systems, that is, to the 
differential-algebraic systems. 



2. Second field equation 

Similarly, also for the second field equation (119), we find only two independent components. Both vanish by the 
assumption of vanishing matter spin r a p. Thus, 

K {2w 6 R'(t) +n 3 X'(t) + [(J, 3 v(t) + Aw 6 u(t)} R{t) + 2 [p 3 u(t) - w 3 v(t)] X{t)} 

- I2g[2m + u{t) + m x v(t)] =24KQe 4e 2 j{e 3 jT 01 )} =0, (170) 

k {fi 3 R'(t) -2w 3 X'(t) + 2\}x 3 u{t) - w e v{t)]R(t) - [Aw 3 u{t) + fj, 3 v(t)]X(t)} 

- 12g [2m x u(t) - m~v(t)] = 24k(? e j[e 2 j(e 3 jt 23 )] =0. (171) 

These are two ordinary linear differential equations (ODEs) of first order in the curvature components R(t) and X(t). 
By suitable linear combinations, we can uncouple the first derivatives of these equations. We add 2w 3 x Eq. (170) to 
fi 3 x Eq. (171) and find, provided A = p\ + Aw 3 w e ^ 0, see (100), 

R'(t) + 2 (u(t) + ^HfcpUw) R(t) - M + H,( t)X(t ) 
Yloi 'k 

^— [2(2m + w 3 + m x p 3 )u(t) + (-m^ p 3 + 2m x w 3 ) v{t)] =0. (172) 

Similarly, we add p 3 x Eq. (170) to — 2wqX Eq. (171) and find, provided jiif. + 4w 3 wq ^ 0, 
X'(t) + 2 (u(t) Ma" 3 ^ 6 ^*)) X(t) + &±Mv{t)R{t) 

_ 12^/k [(2m + n 3 -4m x w e )u(t) + (2m'w 6 + m x /j, 3 ) v(t)] =0. (173) 

Let us introduce, by also using (66), the following abbreviations 

2 Ag 2m~w 6 + m x u 3 2 Ag2m + w 3 + m x u 3 
- :=- S , S , (174) 



Ag m fi 3 — 2m x w 3 2 2m + p 3 — Am x we 

' ~ A 



W 3 ; = — X • (1 ?5 ) 



The constants uJq and uj\ are of even parity, whereas uj\ and w| are of odd parity. This allows us to give the components 
of the second field equation (172) and (173) a more compact and transparent form, 

R\t) = 6w 2 u(i) - 3u%v(t) - 2u(t)R(t) + ^ [(/x§ + Awf) X(t) - 2p 3 {w 3 - w 6 )R(t)] , (176) 
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X'(t) = 3uju(t) + 3w%v(t) - 2u(t)X(t) 



A 



[((4 + Awl) R{t) - 2^3(^3 - w 6 )X(t)] 



(177) 



The choice of signs in (174) and (175) will be motivated in the next subsection, for the moment they are just a 
short-hand notation for certain non linear combinations of the fundamental coupling constants of the theory. 
It may be a bit more transparent, to put (176) and (177) into a matrix form and to reinsert A: 



d (R(t)\ 
dt \X(t)J 



v(t) 



- 2 



u(t) + fl 3 



w 3 ~ w 6 



u(t) 
v(t) 

l4 + *>4\ (R(t)\ 
/i| + 4w 3 w 6 \p 2 3 + 4w| J \X(t)J 



p,% + Aw 3 w 6 



v(t) 



R(t) 
X{t) 



(178) 



We have eight variables but only seven relations between 
them. However, we still have to choose an appropriate 
equation of state p = p{p). In cosmological models for 
late times, p(t) ps is a widespread assumption. 



E. Rearranging the field equations into first order 

form 

For certain purposes (in particular numerical evolu- 
tion) it is more convenient to replace the 3 dynam- 
ical second order equations for the gauge potentials 
by 6 first order equations for the observable quantities 
{a, H, u, v, R, X}. From (139,152,162) one can obtain the 
first order set 

a'(t) = H(t)a(t), (179) 
H'(t) = -2H 2 {t)-J^ )+ JL{«[>( t )-3p(t)] 

+4A - (a - ^) R(t) - (b + a 2 )X(t) 

+ ^(a 2 -4a 3 )v 2 (t)}, (180) 

u'(t) = u 2 (t) -3H(t)u(t) + —\K[p(t)-3p(t)] 

3a 2 I- 

+4A - a Q R(t) - (b + <J 2 )X(t) 
-3a 3 « 2 W} , (181) 

v'(t) = - l -X(t)-v(t)[iH(t)-2u(t)]. / (182) 

along with Eqs. (176, 177) for R'(t),X'(t). In addition 
to the dynamical geometric variables, these equations 
also include the material energy density and the pressure, 
p(t),p(t). The material energy density p(t) is related to 
the dynamical variables by (167) — a relation which could 
be used to eliminate it from the system. The energy den- 
sity and pressure are necessarily related by 

p'(t) = -3H(t)[p(t)+p(t)}. (183) 

This relation follows, on the one hand from the ba- 
sic Noether symmetry conservation law applied to the 
source, and, on the other hand can be derived directly 



from (167) using the system of 6 first order equations 
(176, 177, 179-182). 

The above system needs to be supplemented by an ap- 
propriate equation specifying a relation for p(t) in terms 
of suitable dynamical variables. In GR such a relation 
is often taken in the form of a fluid equation of state 
p = p(p). One could also use such an assumption for 
our PG model or, more generally, one could consider any 
specific relation of the form 

p = p(a,H,u,v,R,X) (184) 

to reduce the dynamical equations of the model to a 
closed system of 6 nonlinear coupled first order ordinary 
differential equations (1st order ODEs) describing the dy- 
namics of 3 geometric degrees of freedom. Alternatively, 
if one had a relation of the form 

p'(t) = f\p(t),p(t),u(t),v(t),---}, (185) 

this would be sufficient for integrating the system (176, 
177,179-183,185). 

To investigate more general models in the presence of 
torsion, one cannot just prescribe a simple equation of 
state. One could consider some explicit source fields and 
their dynamics; these sources would determine p(t),p(t). 
Also one could relax the assumption of vanishing source 
spin density. In follow up work some systems will be 
presented which might be more stable numerically. 

The prescription of an (algebraic) equation of state re- 
duces the phase-space to one of 6 dimensions. Systems of 
ODEs with algebraic constraints are usually called differ- 
ential algebraic equations (DAEs) or also descriptor sys- 
tems 7 . For the numerical evaluation we would have a 
DAE of the form 6© 1, that is, 6 ODEs of first order and 
one algebraic equation which makes the whole system of 
equations determinate. 

Numerical simulations of those systems, including the 
case u(t) ~ v(t), will be discussed in detail in a continua- 
tion of this paper. For the subcase of the Shie-Nester-Yo 



See the link http://www4.ncsu.edU/eos/users/s/slc/www/ 
RESDESCRIPT/resdescript.html: "Usually the term DAE refers 
to systems of ordinary differential equations F(x',x,t) = with 
the Jacobian of F with respect to x' being singular." 



21 



Lagrangian (53), such computations have been already 
done by Li, Sun, and Xi [72-74]. 



F. Acceptable choices of signs 

Concerning the acceptable choices of signs for the pa- 
rameters: at the Lagrangian level of analysis we know of 
only one necessary requirement, namely, in order to sat- 
isfy the principle of least action it is absolutely necessary 
to take kinetic energy terms — here meaning specifically 
the quadratic-in-time-derivative terms — to have a posi- 
tive coefficient. 

For our model this is sufficient to fix the signs for the 
quadratic-in-curvature terms, since it turns out that the 
scalar curvature and pseudoscalar curvatures are each 
linear in the time derivatives of certain connection co- 
efficients. Hence physically in Eqs. (105) one must take 
only the case Ai < 0, A2 < 0. 

Regarding the quadratic torsion terms, the situation 
is not so simple. For our cosmological model at least, 
from (139), (143), and (145), it can be seen that only to 
contains a time derivative of a gauge potential (specifi- 
cally, R'). Thus, by the least-action requirement on this 
quadratic-in-time-derivative kinetic term one must take, 
from (64), the coefficient < 0. Consequently, since by 
convention ao > 0, one should require m + > 0. 

Beyond these considerations, one can ascertain which 
constraints are physically appropriate for the parameters 
only from a detailed analysis of the equations of motion. 
This is left to future work. 



VI. DISCUSSION AND CONCLUSION 

In this work we introduce systematically the notion of 
even and odd parity terms for the construction of gravi- 
tational field Lagrangians in the context of the Poincare 
gauge field theory (PC). Exploiting the theory of alge- 
braic invariants, a Lagrangian results that is at most 
quadratic in the field strengths torsion and curvature. 
Here, we rigorously include interaction terms of even and 
odd parity form, like those of (vector torsion x axial 
vector torsion), that is, V A *A, and (scalar curvature 
x pseudoscalar curvature), that is, R A *X. To obtain 
some insight into the dynamics of those 'shadows' in a 
physically realistic model, we constrained ourselves to a 
model containing only scalar and vector parts and their 
corresponding axial versions. The Vbhn Lagrangian (64) 
can be viewed as a generalization of the recently pre- 
sented Vsny' Lagrangian (54) of Chen et al. [1']. In 
light of the difficulties caused by non-linearities [c -40], 
our model (aside from adding a couple of unimportant 
terms quadratic in the non-dynamic WT a torsion com- 
ponents) may well be the most general PG model that 
can be expected to have a dynamical connection with 
well behaved dynamics. 



From a theoretical point of view, the inclusion of ad- 
ditional interaction terms of odd parity character could 
explain some empirical facts we are faced with in cosmol- 
ogy. Besides the usual handling of even parity functionals 
of the field variables, we treat those of odd parity charac- 
ter on the same footing. This may open the discussion to 
explain the empirical imbalance of matter and antimatter 
on a cosmological scale and other related questions that 
are still open. 

Empirically, the inclusion of additional parameters be- 
yond those of the model (54) will enhance the capacity 
to account for the accelerated universe observations in 
terms of dynamical geometry — dark energy could be a 
PG dynamical connection. It is noteworthy that with 
the new pseudoscalar cross coupling parameters the ac- 
celeration of the universe can be more directly influenced 
by the 0~ mode, see (180), which is known to also couple 
to fermion spin. 

In this work, we present for the first time (as far as we 
know) the notion of the diagonalization of a Lagrangian. 
This identifies certain special parameter combinations 
of the primary coupling constants that are expected to 
play important roles in future studies of the dynam- 
ics of our model, and leads to the recognition of cer- 
tain conditions on the set of primary coupling constants 
{<Z2, 03, W3, u>6> M3} such that the (diagonalized) kinetic 
energy matrix T has strictly positive entries. A working 
hypothesis is that these conditions are needed to have a 
well-defined propagation of massive modes. 

According to the diagonalization, the irreducible pieces 
( 2 *>T a and ^T a can be associated with the two four- 
vectors and a M of even and odd parity character. For 
these vectors to have proper evolution in time, it is ex- 
pected that certain signature properties (conjectured to 
be positive) for the corresponding eigenvalues of the ki- 
netic energy matrix are necessary. The irreducible cur- 
vature pieces i ^>R a P and ( 6 )R a P are essentially a scalar 
R and a pseudoscalar X, respectively. The proper choice 
of parameters is such that the associated kinetic matrix 
has positive eigenvalues. 

In our model, we noted that in the general PG weak 
gravity sector, mediated by the coframe $ Q , the associ- 
ated field strength, the torsion, could carry modes of spin 
2, of spin 1, and of spin (each with even and odd par- 
ity). Restricting to the even parity terms, this is similar 
to Bekenstein's TeVeS (tensor- vector-scalar) theory [97]. 
However, in our case the different modes are carried by 
the torsion alone. There is no need for any other scalar 
or vector fields. 

For strong gravity, mediated by the Lorentz con- 
nection r Q/3 , in our model we found spin of both 
parities. This restriction to zero spin modes in due to 
our simple Lagrangian (64) in which only the scalar 
and pseudo-scalar pieces of the curvature were allowed. 
Straightforward generalizations are possible. In metric- 
affine gravity (MAG), even the inclusion of spin 3 modes, 
see [17], is possible — and all of this on the basis of a 
Riemann-Cartan or metric-affine geometry of spacetime, 
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respectively. 
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